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ABSTRACT 

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude ferm/-LAT 
sources, we have discovered the middle-aged pulsar J2030H-3641, associated with IFGL J2030.0H-3641 
(2FGL J2030.0H-3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have 
obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the 
entire Fermi mission. With a rotation period of 0.2 s, spin-down luminosity of 3 x lO^'* erg s"', and character- 
istic age of 0.5 Myr, PSR J2030H-3641 is a middle-aged neutron star with spin parameters similar to those of 
the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, 
primarily because of its much larger distance, as suggested by the large integrated column density of free elec- 
trons, DM = 246 pc cm"^. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, 
to four geometrical models of magnetospheric emission, and while none of the fits have high significance 
some are encouraging and suggest that further refinements of these models may be worthwhile. We argue 
that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsi- 
ble for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive — 
PSR J2030H-3641 would have been found blindly in gamma rays if only > 0.8 GeV photons had been consid- 
ered, owing to its relatively flat spectrum and location in a region of high soft background. 
Subject headings: gamma rays: observations — pulsars: individual (PSR J2030H-3641) 



1. EVITRODUCTION 

Three years after the launch of the Fermi Gamma-ray Space 
Telescope, pulsations have been detected with its Large Area 
Telescope (LAT) from at least 100 rotation-powered neutron 
stars (e.g., Abdo et al. 2010d), increasing by more than an 
order of magnitude the number of pulsars known to emit 
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above 0. 1 GeVp] Nearly 20 of these are milUsecond pulsars 
(MSPs) known independent from Fermi, for which rotational 
ephemerides obtained from radio observations were used to 
fold th e sparse gamma-ray photons (e.g., |Abdo et al.|2009a[ 
|2010a| l. Another 34 have been discover ed in direct pulsa- 
tion searches of the gamma-ray data (e.g.,|Abdo et al.|2009b[ 
IPletsch et al.pOTT]|SazParkinson et al.|2010| l. 

Among the known non-MSP gamma-ray pulsar sample, 
about 75% have large values of spin-down luminosity (E > 
10^''' erg s"'); the remaining 17 have lower E and are older, 
with characteristic age = P/(2P) ~ 0. 1-1 Myr. Five of these 
are also radio emitters (two are exceedingly faint an d were 
detected onl y after th e gamma-ray discoveries; jCamilo et al.| 
2009b; Plet sch et al.|2 01 1 ). Almost all of these middle-aged 
pulsars are also known or expected to be nearby, at < 1 kpc, 
although obviously many more similar pulsars exist at larger 
distances. In principle several of them could be detected in 
gamma rays despite a relatively low E, because the efficiency 
for conversion of rotational kinetic energy into gamma rays, 
7] = L^/E, apparently increases with dec reasing E and may 
approach 100% for E w 10^^"^'* erg s"'(see |Abdo etal.|2010d 
[Srons 1996; Muslimov & Harding 2003). 

In order to improve our understanding of pulsar emis- 
sion mechanisms and evolution, it is important to iden- 
tify more gamma-ray-emitting neutron stars, particularly in 
under-represented groups such as middle-aged pulsars. The 
Fermi LAT First Source Catalog (IFGL), based on 1 1 months 
of surve y data and containin g more than 600 unidentified 
sources (|Abdo et al.||2010b|l, pr ovides a path towa rd such 

m 



discoveries (as does now 2FGL; Abdo et aLl|2011a] l. Many 



" Seehttps://confluence.slac.stanford.edu/display/GLAMCOG/Public+List+of+LAT- 
Detected+Gamma-Ray+Pulsars for an up-to-date list. 
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unidentified LAT sources are being surveyed within the con- 
text of a "pulsar search consortium" that aims to efficiently or- 
ganize this c ollaborative radio and gamma-ray work (see, e.g., 
[Ray & iS kinson 2011). Spectacularly, many radio MSPs 
have recently been discover ed in suc h searches at high Galac- 
tic latitudes (\b\ > 5°; e.g ., |Cognard et al.,,2011^ ,Keith et al., 
2011 Ransom et al.|20lT i. Here we report on the hrst pulsar 
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to be discovered at low Galactic latitude that is responsible 
for a formerly unidentified IFGL source. PSR J2030H-3641 
in IFGL J2030.0H-3641 (also 2FGL J2030.0H-3640), discov- 
ered with the NRAO Green Bank Telescope (GBT), is also the 
first non-MSP gamma-ray pulsar identified in radio searches 
of IFGL sources. 

2. OBSERVATIONS AND RESULTS 
2.1. Radio Searches 

Consistent with the properties of known gamma-ray pul- 
sars, our first IFGL search targets are non-variable and have 
spectra consistent with exponentially-cutoff power laws. In 
this initial small radio survey at the GBT we aimed to search 
along the Galactic plane, indicating a relatively high search 
frequency (> 1 GHz) in order to minimize deleterious propa- 
gation effects in the interstellar medium and the background 
temperature due to Galactic synchrotron emission. At the 
2 GHz frequency used, the GBT beam is 3' HWHM, and the 
IFGL sources to be searched should therefore, ideally, have 
(95% confidence level error radii) below that, after ac- 
counting for the combined statistical and estimated system- 
atic positional uncertainties. These are strict criteria, and we 
searched a total of only three sources. For each we provide 
(a,S,Ag,ri)s), which are respectively the R.A. and decl. of 
the observed position, its offset from the IFGL position (we 
observed in late 2009, based on an interim version of the cata- 
log), and rgs for the gamma-ray source: IFGL J0224.0h-6201c 
(36.°00,62?04, l',4'), where the "c" indicates that the mea- 
sured properties of this source (including position) may not 
be reliable; IFGL J1746.7-3233 (266.°70,-32?61,3',3'); and 
IFGL J2030.0H-3641 (307P52, 36P68, 1', 3'). 

We observed each source for 1 hr using the GUPPI spec- 
tromete:[^ at a central frequency of 2 GHz. Each of 512 
polarization-summed frequency channels, spanning a band- 
width of 800 MHz, were sampled every 0. 164 ms before writ- 
ing to disk. The flux density limit of these searches, converted 
to a standard search frequency of 1 .4 GHz, was approximately 
0.03 mJy, for an assumed pulsar duty cycle of 10% and spin 
period larger than a few tens of milliseconds. We analyzed 
the data with standard pulsar search techniques implemented 
in PRESTO (Ra nsom 2001), selecting trial dispersion mea- 
sures so as to mamtam ideal sensitivity up to twice th e maxi- 
mum Gal actic DM predicted in each direction by the ICordes 
|& Lazio| (jS)02) electron distribution model. The data sets 
retained significant sensitivity to (possibly binary) MSPs up 
to DM ss lOOpccm"^. We therefore did the analysis in two 
passes, including one where we allowed for a modest amount 
of unknown constant acceleration (parameterized in PRESTO 
by zmax = 50). In the IFGL J2030.0H-3641 data, collected on 
2009 November 27, we discovered on Christmas day an obvi- 
ous pulsar with period P = 0.200 s and DM = 247 pc cm"^. 

2.2. Radio Timing and Polarimetry of PSR J2030+3641 



' https://wikio.nrao.edu/bin/view/CICADA/GUPPiUsersGuide 
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Fig. 1.— Polarimetric observation ofPSRJ2030+3641 at 1.5 GHz with tlie 
GBT. In the lower panel the total-intensity profile is represented by the black 
trace, red corresponds to linear polarization, and blue to circular. The profiles, 
based on a 4.5 hr integration, are displayed with a resolution of 256 bins and 
arbitrary phase. In the upper panel the position angle of linear polarization 
is displayed for bins in which the linear signal-to-noise ratio exceeds 2.2, 
and the PAs have been rotated to the pulsar frame using the measured RM = 
+514radm"^. The full pulse phase corresponds to P = 0.200s. 

We began regular timing observations of the new pulsar at 
GBT on 2010 January 8. Observing parameters were identical 
to those used in the search observation, although with much 
shorter integration times, typically 5 minutes. In this man- 
ner we measured pulse times of arrival (TOAs) on 28 sepa- 
rate days through 201 1 April 11. Using these with TEMPC[^ 
we obtain a phase-connected rotational ephemeris with 0.3 ms 
rms residual. After 5 months of radio timing, the solution was 
sufficient to yield an initial detection of gamma-ray pulsations 
(see Section [273] i. The DM measurement was improved with 
the aid of one additional observation centered at 1 .5 GHz, and 
implies an uncertainty in the alignment of radio and gamma- 
ray pulses of only 0.7 ms. 

In order to determine the polarization characteristics of 
PSR J2030H-3641, we have used GUPPI to make three flux- 
calibrated full-Stokes observations of the pulsar, one each at 
central frequencies of 0.82 GHz (lasting for 2.3 hr), 1.5 GHz 
(4.5 hr), and 2GHz (l.Ohr). W e analyzed the data with 
PSRCHIVE (|Hotan et an^004|, and show the resulting 
1 .5 GHz profiles in Figure llj ITie profiles at the other two 
frequencies are comparable, with a possible slight hint of 
interstellar scattering visible at 0.8 GHz. The Faraday rota- 
tion measure is RM = +514 rad m"^, which implies an average 
electron-weighted Galactic magnetic field strength along the 
line of sight of 2.6 /iG. The radio spectral index, based on 
the three flux density measurements, is a = -1 .7 ± 0. 1 (where 
S^ (X v"; see Table [ij. 

2.3. Gamma-ray Study of PSR J2030+3641 

We used "Pass 6 diffuse"-class Fermi LAT events (having 
the highest probability of being gamma-ray photons), exclud- 
ing those with zenith angles > 100° to minimize gamma-ray 
emission from the Earth's limb. The events were analyzed 
using the LAT Science Tool^^ and photon phases were cal- 
culated using TEMP02 ( |Hobbs et al.|2006| l with the f ermi 
plugin ( ,Ray et al.|201 1) . 

" http://www.atnf.csiro.au/research/pulsar/tempo 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 
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Fig. 2. — Phase-aligned radio and gamma-ray profiles of PSR J2030+3641. 
Bottom panel: discovery observation at the GI3T. Middle panel: Fermi LAT 
0.5-0.8 GeV profile (significant LAT pulsations are not detected below this 
energy band). Top panel: > 0.8 GeV and > 3 GeV (lower) LAT profiles. 
Each pulse profile is repeated in phase, and displayed with 30 bins per period 
in gamma rays and 128 bins in radio. The dashed line gives the level of 
unpulsed emission in the aperture as estimated from the light curve fits. 



We first folded photons detected from IFGL J2030.0H-3641 
based on the initial radio timing solution for PSR J2030H-3641 
(Section |2.2|i. Selecting photons using guesses for the mini- 



mum energy and radius of interest for extraction, we made 
a clear detection of pulsations. Optimal cuts in this re- 
gion of very high background were subsequently found with 
E > 0.8 GeV and within a 0?5 radius. We then obtained 20 
gamma-ray TOAs from a total of 1402 photons detected be- 
tween 2008 August 4 and 2011 April 20 (49 days of integra- 
tion per TOA), and used TEMP02 to determine a timing solu- 
tion with an rms of 4.0 ms. Despite the much higher precision 
of the individual radio TOAs, the best overall ephemeris is 
obtained by combining the gamma-ray and radio TOAs, ow- 
ing to the longer span of the former. The parameters of this 
timing solution, in which a small amount of "timing noise" is 
parameterized by i/ (where iy= 1 /P), are presented in Table[T] 
Figure l2] shows the gamma-ray pulse profiles in three en- 
ergy bands^ along with the radio discovery profile. In order to 
optimize each pulsed signal-to-noise ratio, we selected pho- 
tons within 0?8, 0?5, and 0.°3 of the timing position, respec- 
tively for the 0.5-0.8 GeV, > 0.8 GeV, and > 3 GeV bands. 
For the middle energy band, the bin-independent H-test (jde] 
Jager & Busching 2010), with a value of 270, indicates enor- 
mous significance. For each of the other two bands, the H- 
test values of approximately 40 correspond to pulsation sig- 
nificances over 5(7. The light curve at > 0.8 GeV consists of 
two (conceivably three) overlapping components as observed, 
e.g for PSRs J1709- 4429, J1057-5226, and the CTA 1 pul- 
sar ( Abdo et al. 2010d), which have peak separation A w 0.25. 
While the statistics are poor, the > 3 GeV profile is consistent 



with such a description. At 0.5-0.8 GeV, only one peak is 
visible, which corresponds in phase with the first peak visi- 
ble at higher energies. No pulsed signal was detected below 
0.5 GeV, despite searches using different apertures. We fit 
the unbinned gamma-ray data above 0.8 GeV with two two- 
sided Gaussians, to account for the steep leading and trail- 
ing edges of the pulsed emission. The first peak is offset 
in phase from the main radio peak hy 6 = 0.26 ± 0.02 with 
FWHM = 0.11 ±0.04, while the second peak is offset from 
the first by A = 0.30!|J:!]^, with FWHM = 0.18 ± 0.05. We 
have included systematic error from the modest uncertainty 
in DM and our choice of model. 

2.3.1. Spectrum 

To characterize the spectrum of PSR J2030H-3641, we an- 
alyzed Pass 6 diffuse events collected between 2008 August 
4 and 2011 April 20 selected to have reconstructed energies 
between 178 MeV and 100 GeV and reconstructed positions 
within 10° of the timing position of PSR J2030H-3641. We 
excluded the 100-178 MeV energy bin because the rapidly 
increasing effectiv e area in this bin exa cerbates the effects of 
energy dispersion ( Atwood et al.|2009 1. We removed events 
with a reconstructed zenith angle > 100 and also filtered data 
in which the aperture approached the limb of the Earth; these 
cuts remove most of the contribution of the bright gamma-ray 
emission produced by cosmic rays interacting in the Earth's 
atmosphere. To further enhance the signal-to-noise ratio of 
PSR J2030H-3641, we selected events with rotation phase 
0.22 < <j) < 0.64 (see Figure l2|. We verified that there was 
negligible emission outside ortnis phase range, so the best-fit 
values may be taken to represent the phase-averaged magne- 
tospheric emission. 

PSR J2030H-3641 lies in the source-crowded Cygnus re- 
gion, and extracting its spectrum requires careful model- 
ing of both neighboring point sources (including the three 
bright gamma-ray pulsars J2021H-3651, J2021H-4026, and 
J2032H-4127) and the diffuse background. We included all 
point sources in a p reliminary version of the 2FGL catalog 
(Abdo et al. 201 la[ ) within 15° of the pulsar, as the broad 
point- spread function (PSF) at low energy allows sources 
away from the aperture boundary to contribute to the observed 
counts. Identified pulsars were modeled with a power law 
with exponential cutoff, dN/dE cx (£'/GeV)"'"exp(-£'/£'c)- 
In the fits reported below, we varied all parameters for sources 
within 7° of the aperture, while sources outside of the aperture 
were held fixed at the catalog values. We modeled the Galac- 
tic and isotropic diffuse backgrounds using, respectively, 
the gll_iem_v02_P6_Vll_DIFFUSE diffuse map and the 
isotropic _iem_v02_P6_Vl 1 _DIFFUSE tabulation employed 
for the IFGL spectral analysis; these have been rescaled from 
those constructed with the P6_V3_DIFFUSE instrument re- 
sponse functions (IRFs). We allowed the normalizations of 
the diffuse models to vary in the fit. 

We perfo rmed spectral fits with the pointlike application 
( Kerr|20lT l. This software bins the events in both energy and 
position, using HEALPi^jp*] to provide position bins whose 
energy-dependent extent is always small compared to the PSF, 
providing a compact representation of the data with minimal 
information loss. We verified the results with the Science- 
Tool gtlike. The best-fit values are reported in Table [T] and 
include a power-law photon index F = 1 . 1 and cut-off energy 

^' http://healpix.jpl.nasa.gov/ 
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Ec = 2.0 GeV. These values apply to the profile as a whole, 
although Figure |2] suggests that the second gamma-ray peak 
has a harder spectrum than th e first, as is ob served for many 
gamma-ray pulsars (e.g., Abdo et al.|201 lb| l. 

In addition to the rescaled IFGL diffuse model, we also per- 
formed the fit with an improved model internal to the LAT col- 
laboration and with "Pass 7" datij^using the P7SOURCE_V6 
IRFs and the appropriate diffuse model (accompanying the 
2FGL catalog). These values were consistent with those re- 
ported in Table [T] and the scatter has been used to estimate 
systematic errors on the parameters. We have also included 
uncertainty in the instrument's effective area in the system- 



atic err ors through the use of "bracketing" IRFs ( |Abdo et al. 
|2009c| l. 

2.4. X-ray Observations 

A short Swift ( [Gehrels et al.||2004| observation was ob- 
tained on 2010 May 1 1 to search for an X-ray counterpart to 
PSR J2030H-3641. No emission is o bserved at the puls ar po- 
sition in this 3.8 ks X-ray Telescope ( Burrows et al.|200 5 ) PC 
mode observation, with a 3a upper limit of 1.7 x I0~^ cts s"' 
(0.5-8 keV). Assuming a power-law spectrum with p hoton 
index F = 1.5 for a pos sible pulsar wind nebula (e.g., Kar- 
Igaltsev & Pavlov] |2008| l absorbed by a column with Nh = 
7.4 X 10^' cm""^ (assuming one free electron for every 10 neu- 
tral hydrogen atoms along the line of sight), this corresponds 
to an unabsorbed flux limit of fx < 1 .3 x lO"'-' erg cm"^ s"', or 
luminosity Lx < 1.6 x lO-'H'^/l kpc)^ erg s"'. For a distance 
of 2 kpc (see Table[l|, this corresponds to Lx/E <2x 10~^, 
which is a reasonably stringent limit, although it would not be 
a surprise if the actual X-ray flux of this pulsar were about an 
order of ma gnitude below the current limit (see, e.g., |Marelli| 
let al.|201I| i. 

3. DISCUSSION 

The gamma-ray pulsations from PSR J2030H-3641 unam- 
biguously identify it as the origin of the gamma rays from 
IFGL J2030.0H-3641. We note also that a candidate TeV 
source possibly dete cted by the M ilagro Gamma-Ray Ob- 
servatory ("C2" in Abdo et ar||200 7) is positionally coinci- 
dent with PSR J2030H-3641. Its faintness and location near a 
very bright source precludes fu rther investig ation based only 
on published results (see also |Abdo et~ar |2009d), but the 
relatively small E and relatively large age and distance of 
PSR J2030H-3641 make an association unlikely. With a DM- 
derived distance of 8 kpc, PSR J2030H-3641 is located toward 
the Cygnus region, where until recently very few pulsars were 
known and where the pulsar distance scal e is highly uncer- 
tain. Recent discoveries suggest that the Cordes & Lazio 
( |2002 1 model greatly overpredicts distances in this direc- 
tion: PSR J2032H-4127, located 5° from PSR J2030H-3641, 
with a DM-derived distance of 3.6 kpc, is thought to be 



located at half that distance ( |Camilo et al. |20(j9b| l; and 
PSR J2021H-3651, less than 2° from PSR J2030H-3641, with 
a model distance of 12 kpc, is more likely located at 2^ kpc 
( |Abdo et al.||2009el l. Also, the RM for PSR J2030H-3641 is 
large and positive, very similar to that of PSR J2021H-3651, 
and this is suggestive of a location closer than the Perseus arm, 
at 6 kpc in this direction (see discussion in Han et al. 2006). 
Overall, we judge that a hkely distance for PSR J2030H-3641 
is in the range 2-4 kpc. 



3.1. Geometrical Constraints from Radio Profile 

[Johnston & Weis berg| ( |2006[ ) noted that the profiles of pul- 
sars with high E share many common features. In particular, 
the profiles with two components generally show (i) two com- 
ponents with equal width, (ii) the trailing component stronger 
than the leading component, (iii) close to 100% linear polar- 
ization, (iv) significant circular polarization under the trail- 
ing component, and (v) flat position angle swings which ap- 
pear to steepen at the far trailing edge of the profile. The 
profile of PSR J2030H-3641 conforms in some ways to these 
generalizations. The profile is symmetrical with the trailing 
edge brighter The PA swing is remarkably flat with a hint of 
steepening at the trailing edge. However, although the leading 
component is highly polarized, the tr ailing edge is not, unlike 



the ma jority of the high E pulsars. Weltevrede & Johnston 
( |2008p pointed out that the transition between pulsars with 
low linear polarization and those which are highly polarized 
occurs at £ - lO^'*-'' erg s"' . PSR J2030H-3641 could be an ex- 
ample of a pulsar in transition from high to low polarization. 



The interpretation of these profiles by Johnston & Weisberg 
( |2006 1 was that the magnetic pole crossing occurs at the (sym- 
metry) center of the profile. The PA swi ng is significantly 
shifte d to later phase because of aberration (Blaskiewic z et al.| 
I99I| l, leading to emission heights close to 1000km. This 



may also be the case in PSR J2030H-3641. If we assume that 
the inflexion point of the PA swing is the rise in the PA curve 
some 20° after the profile center in Figure [T] this implies an 
emission height of ^ 800 km (close to 10% of the light cylin- 
der), similar to those seen in other high E pulsars. We also 
note that a height of ^ 800 km would in turn imply an overall 
pulse width of some 50°, similar to the observed value. This 
further suggests that the inclination between the magnetic and 
rotation axes cannot be too far from orthogonal. 

3.2. Light Curve Fitting and Model Implications 

We can use the gamma-ray light curve and polarimetric 
measurements to test the predictions of magnetosphere beam- 
ing models. Such tests are most constraining when the mag- 
netic inclination angle a and the viewing angle ( are well 
determined. The polarization data provide some geometri- 
cal constraints. Since the lower frequency data may be af- 
fected by interstellar scattering, our primary constraints are 
based on the 1 .5 GHz data. We use PA values from the na- 
tive 512 observational phase bins with ctpa < 10° and phases 
-0.1 < (p < 0.05 (relative to the intensity maximum; see Fig- 
ure [T]i, and fit to a standard rotating vector model (RVM) po- 
larization sweep. The maximum sweep rate is d^'/d0|max 
= 1 .4 and one obtains good fits in a band of the a-C plane, 
with a best fit = 37.5 for 55 degrees of freedom (suggesting 
the PA errors are somewhat overestimated). Figure |3] shows 
the RVM fit confidence regions in this p lane, after correcting 
for th e RVM sign convention problem ( jEverett & Weisberg] 
|2001| l. The phase of maximum sweep is not well constrained, 
but the best fits place it well after the radio peak at </> w 0.06. 
The fit implies a < 140°, with best fit values < 90°. If we si- 
multaneously fit the 1.5 GHz points and the -0.1 < cj) < 0.02 
0.8 GHz data, larger a > 90° are preferred, but this preference 
is controlled largely by the trailing edge of the sweep where 
scattering di stortions may be present. 

Following Romani & Walters ( 20101), we next compare the 
gamma-ray light curve to three geometrical models of outer 



See http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitooIs/ 



magnetosphere emission. The first, "Outer Gap" (OG; Cheng 
[eral.,1986) , model posits emission from vacuum dipole field 
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Fig. 3. — Surfaces of goodness of fit (xs) to the E > 0.8 GeV light curve for the OG (left) and SL (right) outer magnetosphere beaming models. Dai'k colors 
indicate lower (better) values of the fit statistic. The Ax^ = +l,+2,+3 contours for the RVM fit to the 1.5 GHz polarization data are shown in green. 
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Fig. 4. — Model light curves (red) at the best-fit a and ( consistent with 
the RVM constraints compared to the E > 0.8 GeV LAT light curve (black 
histograms). The offset tp indicates the phase of the closest approach of the 
magnetic dipole axis at the star surface. The SL curves have been offset by 
+8() counts/bin. 

lines radiating from the null charge surface to the light cylin- 
der at radius /?lc- In this model, the heuristic efficiency scal- 
ing r] w (10^^ergs-V£)i/2 is assumed. For PSR J2030H-3641, 
the value 77 = 0.18 was used for the gap width w in the OG 
calculations. The second model, "Two Pole Caustic" (TPC; 



|Dyks & Rud ak 2003), has emission in the vacuum open zone 
from the neutron star surface out to 0.75 /?lc from the rota- 
tion axis or an altitude of 7?lc, whichever is lesser. More re- 
cently, numerical models of force-fre e magnetospheres have 
been computed by Bai & Spitkovsky (2010) and they have ar- 
gued that the pulsar gamma-ray beam can be created in the 
"Separatrix Layer" where the open zone current sheets over- 
lap. The emission in this SL model is assumed to extend to 
1.5 .Rlc- We assess the match to the data by plotting a good- 
ness offit Xi surface in the a-C, plane (see Romani & WaF] 
rs||2010| l. Using the E > 0.8 GeV light curve, the best-fit 
OG models (with best X3 = 3.9 at a = 133°, C = 110°) lie in- 
side the 3(7 1.5 GHz fit contours and at the maximum of the 
combined 1.5 GHz and 0.8 GHz fit. A small second region of 
adequate but poorer fits (xi = 7) is also present within the 1 cr 
RVM contour at smaller a = 63° (C = 60°). The SL model 
has a similar small region of good fits within the 1 cr 1.5 GHz 
RVM contour (^3 = 5 at a = 50°, C = 48°). The TPC model 
predicts a large unpulsed component resulting in much poorer 
fits (not shown), with a best value X3 = 19 in the RVM-allowed 
region; this is at a = 26°, C = 25°, giving a light curve with 
little pulsed flux. 

Figure[4]shows the best-fit light curves in the RVM-allowed 
region for the OG and SL models, compared with the data. 
The present models do not include the full radiation physics, 
employing simple uniform emissivity throughout the active 
zones. However, it appears that both high altitude models can 
produce reasonable matches to the LAT data. With improved 
geometrical constraints it may be possible to distinguish be- 
tween the OG and SL models. For the former, the fit phase of 
the magnetic axis is at « 0.09, after the pulse and somewhat 
later than the polarization sweep. The SL model prefers a 
magnetic axis much closer to the maximum of the radio pulse. 
In both cases, the lag of the sweep maximum from the pulse 
peak implies significant altitude ^ 0.1-0.15/?lc for the radio 
emission. Such altitudes are also required to accommodate 
the relatively wide radio pulse in the open zone. Modeling 
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Fig. 5. — Reduced values (with 25 degrees of freedom) in a-^ space 
from slot gap gamma-ray light curve fits, marginalized over w, r, and Aci), and 
black curves from the RVM fit to radio polarization measurements, showing 
where these constraints intersect the geometrical light curve para meter s (the 
overall best-fit light curve lies outside these regions; see Section [T2) . The 
plotted values of xi^, are those included in the Markov chains via the algo- 
rithm given in Verde et al. (2003 1. 

such finite altitude effects may provide some additional re- 
strictions on the geometry. We also note that the OG and SL 
models predict somewhat different viewing angles C^; if an X- 
ray pulsar wind torus could be measured, the resulting <^ might 
be used to select the preferred model. 

The classic TPC model fails in these fits at least in part be- 
cause it does not have enough emission at high altitudes. In 
order to improve upon this, we also fit the > 0.8 GeV profile 
of PSR J2030-I-3641 with light curves predicte d from a ge- 
ometrical high-altitude "Slot Gap" model (SG; [Muslimov & 
"Hardin g 2 004]l, an evolu tion of the low-altitude SG of Mus- 
limov & Harding' (2003 ). In the SG geometry, the emission 
ongmates in a gap that extends from the neutron star surface 
at the polar cap to high altitudes near the light cylinder. We fit 
for four model parameters: a, (, gap width w (a fraction of the 
open volume radius of the polar cap, rpyc, in the open vo lume 
coordinate system described in |Dyks & Harding]|2004 ), and 
maximum radius of emission r (in units of i^LcX assuming 
uniform emission rate in the corotating frame. Our resolu- 
tion is 1° in a and 0.01 rove in < w < 0.3, and 0.1/?lc 
in 0.7/?LC < < 2.0/?LC (we assumed a maximum cylindrical 
emission radius of 0.95 /?lc)- The simulated light curves are 
fit to the data via a Mai'kov Chain Monte Cario (MCMC) rou- 
tine that explores the para meter space a nd records regions of 
maximum likelihood, as in |Verde et al.| (j2003 ). 

As noted before, we have some independent constraints on 
the emission geometry of PSR J2030-I-3641 from radio polar- 
ization (cf. Figure [TJ- To make our model fits be consistent 
with the RVM fit, we allowed the model light curve to shift 
with respect to the location of the magnetic pole by no more 
than A(p = 0.1 in phase. Figure |5] shows the a-( space with 
the polarization contours already shown in Figure[3]and color 
"contours" of reduced from the MCMC likelihood fits. The 
polarization contours pass through parameters yielding light 
curve fits with 3 < < 5, the best fit falling within the 

contours occurring at (a, (, w, r, A(f>; x^/25) = (134°, 1 14.5°, 
0.0 rove, 1.4/?LC, 0.08; 3.4). The second best fit that falls 
within the polarization contours is at (54°, 52.5°, 0.01 rove, 
1.3/?LC, 0.105; 4.2). We include this fit because while it has 
a larger x^, it falls within a more constrained region of the 



polarization contours. The absolute best fit, with or without 
restricting A(/), is at (170°, 110.5°, O.Orovc, > 1.1/?lc, 0.017; 
1 .4). This is inconsistent with standard radio cone beam emis- 
sion because the values of a and ( are too far apart for lower 
altitude radio emission to be visible. All of these light curves 
are shown in Figure [6] 

Two of the fits include w = 0, meaning the emission in the 
model originates from a single open field line. Physically, this 
means the emission zone width is smaller than the resolution 
of our simulations, but not really zero. All the fits have small 
widths, requiring higher current of radiating particles in the 
gap (either more primaries or high pair multiplicity) in order 
to generate the observed luminosity. It is also likely that the 
emissivity and width of the emission zone change azimuthally 
and with altitud e, unlike the constant width of our simula- 
tions (see Harding & Muslimov 2011). None of the fits are 
able to reproduce welTThe small middle peak of the observed 
light curve; it is not yet clear whether this is a separate feature 
or part of the second peak. Also, the background is slightly 
over-predicted by the SG model, especially for the parameters 
consistent with the polarization curves. 

None of these (OG, SL, SG) gamma-ray light curve fits are 
particularly consistent with our inferences from the radio pro- 
file and polarimetric constraints. For instance, the latter sug- 
gest a nearly orthogonal rotator (Section 3. Ill, for which there 



are no good gamma-ray light curve fits^ However, the po- 
larimetric data provide limited constraints, owing to the small 
range of pulse longitude for which they are obtained, and our 
inferences rely on some assumptions. As for the gamma-ray 
fits, none are perfect, although the OG and SG fits appear 
slightly better than the SL one. Considering their inherent 
limitations, however, these fits are encouraging, and further 
refinements may improve our understanding of the magneto- 
sphere of middle-aged pulsars like PSR J2030-I-3641. 

For each beaming model, one can extrapolate the aver- 
age pulsar flux on the Earth line-of-sight to the isotropic all- 
sky equivalent. For the best-fit angles consistent with polar- 
ization constraints, the corrections for PSR J2030-I-3641 are 
fn = 0.66 (OG), 0.64 (SL), and 0.78 (SG; neglecting polar- 
ization constraints, /si.sG = 0.45). Using the observed flux 
and the 8kpc DM-estimated distance, we find that r^obs = 
And^ffiFj/E = 10ii|/o. This indicates that any distance 
greater that ^ 2.5 kpc is problematic, unless the emission is 
highly beamed. For the OG/SL/SG beaming ^ 0.7, we re- 
quire d < 3 kpc. The preferred distance for the heuristic OG 
efficiency law (77 = 0.18) is ^ 1.3 kpc. Such distance estimates 
of « 2 kpc are consistent with the discussion at the beginning 
of Section|3] and strongly suggest that PSR J2030H-3641 is lo- 
cated at roughly lOx the distance to Geminga (d = 0.16kpc). 
Improved distance estimates for middle-aged, apparently effi- 
cient pulsars have the potential to constrain or falsify gamma- 
ray emission models. 

3.3. Where Have All the ( Ordinary Radio) Pulsars Gone? 

In a radio search with the GBT of three unidentified 
low Galactic latitude LAT sources, we have discovered 
PSR J2030-I-3641, a relatively distant member of the under- 
represented group of gamma-ray pulsars that are middle- 
aged and have small E. A gamma-ray-only pulsar was sub- 
sequently discovered in one of the two remaining sources, 
IFGL J1746.7-3233 ( |Pletsch et al.1[20TT] l. The third source 
could yet be a gamma-ray pulsar, but is unlikely to be a 
substantial radio emitter beamed towards the Earth (assum- 
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Fig. 6. — Slot gap model light curves compared to the data, (a) The best-fit light curve (blue) with parameters lying within the polarization contours. The LAT 
light curve is shown in black. The vertical dot-dashed line shows the model location of the magnetic pole, shifted by Aiji from zero. The horizontal dashed line 
represents the background count level, (b) Same as (a), for the second-best fit parameters (red), (c) Same as (a) and (b) for the absolute best-fit light curve (green) 
found by the MCMC routine, with no restriction on A(/>. Its parameters do not fall within the radio polarization contours. 



ing we covered its location adequately — in the case of 
IFGL J1746.7-3233, the pulsar is 3:8 offset from our GBT 
search position; see Section |2.1| i: for a fiducial distance of 
4kpc, our sensitivity corresponds to a luminosity limit L1.4 = 
SiAd^ = 0.5 mJy kpc^, which is near the low end of the radio 
luminos ity distribution for detected radio pulsars (see |Abdo| 
|et al.|2010c, Camilo et al. 2009a b, which include the two ex- 
ceptions with much lower luminosities). It would seem that 
we were quite successful, with one related radio discovery 
among three LAT sources searched. However, as we argue 
next, we don't expect such continued success. 

Many more Galactic plane unidentified LAT sources have 
been searched for radio pulsars than the three mentioned here. 
This work is being done at the GBT, Parkes, Nan9ay, GMRT, 
and Effelsberg telescopes by members of the Fermi pulsar 
search consortium, and includes searches of all sources that 
have pulsar characteristics at all Galactic latitudes. So far this 
has resulted in the discove ry of mor e than 30 MSPs, all at high 

201 l[|Keithetal. 2011 ; Ransom, 
and PSR J20^0+364l. While a'detail^ popu- 



latitudes (e .g., |Cognard et al. 
etal.|2011 



lation analysis will have to wait for the publication of all the 
completed searches, the lone example of PSR J2030-I-3641 al- 
ready stands out. Why have more radio pulsars not been de- 
tected among the unidentified Galactic plane LAT sources? 

A point related to this question concerns the very deep ra- 
dio searches done of 34 pulsars discovered in blind pulsa- 



tion searches of gamma rays (e.g.,|Abdo et al. 2009b ; Pletsch 
2011 Saz Parkinson et al. |2010ir"only four of them 



et al. 



were 



bund toal so be radio emitters (see [Pletsch et al.|2011 

|Ray et aL]|201 1\ . Three are exceptionally faint and their de 
tectio n was possible only following the gamma-ray discov 
eries ( |Abdo eraL||2010c| [Camilo et all2009bt IPletsch et al. 



2011 1. The remaining one, PSR J2032H-4127 ( [Camilo et aL 
2009b|l, is like PSR J2030H-3641 in that they are both located 



in the Cygnus region and both are relatively bright, being de- 
tectable in 1 minute with the GBT at 2 GHz. They were not 
known prior to the launch of Fermi simply because no sensi- 
tive surveys had (or have yet fully) been done of the northern 
(equatorial) reaches of the Galactic plane, using the high radio 
frequencies desired to unveil the Galactic disk. Conversely, 
the enormo usly successful 1.4 GHz Parkes multibeam survey 
(e.g., Man chester et al.[[2001 1 did search the entire Galactic 
plane at 260" < / < 50" with sensitivity sufficient to detect or- 
dinary (non-MSP) radio pulsars with Li 4 > 0.5 mJy kpc^ out 
tod>2 kpc. A very large fraction of all potentially detectable 
ordinary radio pulsars in that area of the sky that could also be 
plausible gamma-ray sources were therefore discovered be- 
fore Fermi launch, especially when adding those discovered 
in de ep pointed observati ons of known pulsar wind nebulae 
(e.g., 'Cami lo et aL]|2002| l. In time, several of these pulsars 
were detected by LAT (e.g., [Abdo et al.|2010d|l. Subsequent 
to the Parkes multibeam survey an even deeper Parkes survey 
extended Galactic plane coverage out to Z = 60° (Camilo et 
al., in preparation), and th e on-going PALFA Arecibo survey 
is doing so out to I sa 75° ( Cordes et al.|2006 i. 

In our view, the above suggests that the answer to "why 
have more radio pulsars not been detected among the uniden- 
tified Galactic plane LAT sources?" is that the vast majority 
of ordinary radio pulsars accessible to the current generation 
of telescopes and located within a few degrees of the plane at 
/ < 75° that can yield an appreciable gamma-ray flux at the 
Earth were discovered a long time ago (it is possible that in 
rare but important cases very high distant pulsars will have 
their radio pulses scattered beyond practical detectability). 

Many of the LAT sources that remain unidentified along 
the Galactic plane are surely pulsars — but they may not be 
detectable as radio sources, and should be searched anew in 
gamma rays (some of these may be in binaries, and inacces- 
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sible to current blind searches; see, e.g., Corbet et al.|201 1\ . 
PSR J2030+3641 is bright enough that it would have been 
discovered in blind searches of 18 months of LAT data, if 
the photon selection criteria had been adjusted to take into ac- 
count its flat spectrum (e.g., if only photons with E > 0.8 GeV 
had been searched for pulsations), as we confirmed after dis- 
covery. Spectral analysis of a region can improve the signal- 
to-noise ratio by selecting events based o n the probability that 
they come from the source of interest (see Belfiore 2011| l, and 
applying this technique to the Cygnus region also resulted in 
the unbiased detection of pulsations from PSR J2030H-3641. 
For fainter point sources superimposed on the large and un- 
certain diffuse Galactic background, spectral analysis and lo- 
calization are harder (IFGL J2030.0H-3641 was only Q'A from 
the actual pulsar position), but these and other improvements 
in search techniques are already result ing in the discovery 
of many m ore gamma-ray pulsars (see [Pletsch et al.||20lTt 
Saz Parkinson 2011). As a consequence, the ratio of known 



Earth (e.g., Romani et al.|201 1 1. 



gamma-ray-only to gamma-ray-and-radio ordinary pulsars, 
which is currently slightly under 1 .0, should increase substan- 
tially. 

When will we ever learn all that we can (in the Fermi 
era) about the geometry and emission properties of gamma- 
ray and radio pulsar accelerators? Perhaps when full details 
emerge from the continuing collaborative radio and gamma- 
ray observational and modeling work. The future is bright (at 
some wavelength), but requires publication of all searches, 
including radio non-detections, and consideration of the non- 
detection by LAT of radio pulsars with very large E flux at the 



The GBT is operated by the National Radio Astronomy Ob- 
servatory, a facility of the National Science Foundation op- 
erated under cooperative agreement by Associated Universi- 
ties, Inc. The Fermi LAT Collaboration acknowledges gener- 
ous ongoing support from a number of agencies and institutes 
that have supported both the development and the operation 
of the LAT as well as scientific data analysis. These include 
the National Aeronautics and Space Administration (NASA) 
and the Department of Energy in the United States, the Com- 
missariat a I'Energie Atomique and the Centre National de 
la Recherche Scientifique/Institut National de Physique Nu- 
cleaire et de Physique des Particules in France, the Agenzia 
Spaziale Italiana and the Istituto Nazionale di Fisica Nucle- 
are in Italy, the Ministry of Education, Culture, Sports, Sci- 
ence and Technology (MEXT), High Energy Accelerator Re- 
search Organization (KEK) and Japan Aerospace Exploration 
Agency (JAXA) in Japan, and the K. A. Wallenberg Founda- 
tion, the Swedish Research Council and the Swedish National 
Space Board in Sweden. Support for this work was provided 
by NASA through Einstein Postdoctoral Fellowship Award 
Number PFO- 110073 issued by the Chandra X-ray Observa- 
tory Center, which is operated by the Smithsonian Astrophys- 
ical Observatory for and on behalf of NASA under contract 
NAS8-03060. 

Facilities: Fermi (LAT), GBT (GUPPI), Swift (XRT) 



REFERENCES 



Abdo, A. A., et al. 2007, ApJ, 664, L91 

— . 2009a, Science, 325, 848 

— . 2009b, Science, 325, 840 

— . 2009c, ApJ, 696, 1084 

— . 2009d, ApJ, 700, LI 27 

— . 2009e, ApJ, 700, 1059 

— . 2010a, ApJ, 712, 957 

— . 2010b, ApJS, 188,405 

— . 2010c, ApJ, 711,64 

— . 2010d,ApJS, 187,460 

— . 2011a, ApJS, submitted (arXiv/1 108.1435) 

— . 2011b, ApJ, submitted (arXiv/1 107.4151) 

Aions, J. 1996, A&AS, 120, C49 

Atwood, W. B., et al. 2009, ApJ, 697, 1071 

Bai, X., & Spitkovsky, A. 2010, ApJ, 715, 1282 

Belfiore, A. 201 1, in Radio Pulsars: an Astrophysical Key to Unlock the 
Secrets of the Universe, ed. M. Burgay, N. D'Amico, R Esposito, 
A. Pellizzoni, & A. Possenti, Vol. 1357 (AIP), 44 

Blaskiewicz, M., Cordes, J. M., & Wasserman, I. 1991, ApJ, 370, 643 

Burrows, D. N., et al. 2005, Space Sci. Rev., 120, 165 

Camilo, E, Manchester, R. N., Gaensler, B. M., & Lorimer, D. R. 2002, ApJ, 
579, L25 

Camilo, E, Ng, C.-Y., Gaensler, B. M., Ransom, S. M., Chatterjee, S., 

Reynolds, J., & Sarkissian, J. 2009a, ApJ, 703, L55 
Camilo, E, et al. 2009b, ApJ, 705, 1 

Cheng, K. S., Ho, C, & Ruderman, M. 1986, ApJ, 300, 500 
Cognard, I., et al. 201 1, ApJ, 732, 47 

Corbet, R. H. D., et al. 201 1, IFGL J1018.6-5856: a new gamma-ray binary, 
ATEL 3221 

Cordes, J. M., & Lazio, T. J. W. 2002, preprint ^arXiv:astro-ph/0207T56) 

Cordes, J. M., et al. 2006, ApJ, 637, 446 

de Jager, O. C, & Busching, I. 2010, A&A, 517, L9 

Dyks, J., & Harding, A. K. 2004, ApJ, 614, 869 

Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 

Everett. J. E., & Weisberg, J. M. 2001. ApJ, 553, 341 

Gehrels, N., et al. 2004, ApJ, 61 1, 1005 



Han, J. L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van Straten, W. 

2006. ApJ, 642, 868 
Harding, A. K., & Muslimov, A. G. 2011, ApJ, 726, LIO 
Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 
Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302 
Johnston, S., & Weisberg, J. M. 2006, MNRAS, 368, 1856 
Kargaltsev, O., & Pavlov, G. G. 2008, in American Institute of Physics 

Conference Series, Vol. 983, 40 Years of Pulsars: Millisecond Pulsars, 

Magnetars and More, ed. C. Bassa, Z. Wang, A. Cumming, & 

V. M. Kaspi, 171-185 
Keith, M. J., et al. 201 1, MNRAS. 414, 1292 

Kerr, M. 2011, PhD thesis. University of Washington, arXiv/1 101.6072 

Manchester, R. N., et al. 2001, MNRAS, 328, 17 

Marelli, M., De Luca, A., & Caraveo, P A. 201 1, ApJ, 733, 82 

Muslimov, A. G., & Harding, A. K. 2003, ApJ, 588, 430 

— . 2004, ApJ, 606, 1 143 

Pletsch, H. J., et al. 201 1, ApJ, in press (arXiv/1 1 1 1.0523) 
Ransom, S. M. 2001, PhD thesis. Harvard University 
Ransom, S. M., et al. 2011, ApJ, 727, L16 

Ray, P. S., & Parkinson, P. M. S. 201 1, in High-Energy Emission from 
Pulsars and their Systems, ed. D. F. Ton'es & N. Rea, 37 
(arXiv/1007.2183) 

Ray. P S., et al. 2011, ApJS, 194, 17 

Romani, R. W, Kerr, M., Craig, H. A., Johnston, S., Cognard, I., & Smith, 

D. A. 2011, ApJ, 738, 114 
Romani, R. W, & Watters, K. P 2010, ApJ, 714, 810 
Saz Parkinson, P. M. 201 1, in Radio Pulsars: an Astrophysical Key to 

Unlock the Secrets of the Universe, ed. M. Burgay, N. D'Amico, 

R Esposito, A. Pellizzoni, & A. Possenti, Vol. 1357 (AIP), 48 

(arXiv/1 101.3096) 
Saz Parkinson, P M., et al. 2010, ApJ, 725, 571 
Verde, L., et al. 2003, ApJS, 148, 195 
Weltevrede, P, & Johnston, S. 2008, MNRAS, 391, 1210 



PSR J2030+3641: MIDDLE-AGED RADIO AND GAMMA-RAY PULSAR 



9 



TABLE 1 

Measured and Derived Parameters for PSR J2030+3641 



Parameter Value 



Right ascension, R.A. (J2000.0) 20''30""00?261(4) 

Declination, decl. (J2000.0) 36°41'27"15(7) 

Galactic longitude, / (deg) 76.12 

Galactic latitude, (deg) -1.44 

Rotation frequency, v 4.99678736138(4) 

Frequency derivative, v (i^^) -1.62294(5) X 10" 

Frequency second derivative, i) {s"^) 2.2(6) X lO"^** 

Epoch of frequency (MJD) 55400.0 

Dispersion measure, DM (pccm^^) 246.0(7) 

Data span (MJD) 54700-5567 1 

RMS timing residual (ms) 0.4 

Spin-down luminosity, E (ergs"') 3.2 X 10^** 

Characteristic age, Tc (yr) 4.9 X lO' 

Surface dipole magnetic field strength (Gauss) .... 1 .2 X lO'^ 

Flux density at 2 GHz, S2 (mjy) 0.09 

Flux density at 1 .5 GHz, Si 5 (mJy) 0.15 

Flux density at 0.8 GHz, So.g (mJy) 0.40 

Rotation measure, RM(radm-2) +514±1 

Radio-gamma-ray profile offset, 5 (P) 0.26 it 0.02 

Gamma-ray profile peak-to-peak separation, A (P) 0.30^5J [Jg 

Gamma-ray (> 0. 1 GeV) photon index, V 1 . 1 ± 0.2;p 

Gamma-ray cut-off energy, Ec (GeV) 2.0 it O.S^jj^ 

Photon flux (> 0.1 GeV) (lO^^^cm^^ s"') 3.6it 0.5 ±0.8 

Energy flux (> 0.1 GeV), P^ (10"" ergcm^^s-') . 4.2 it 0.3 it 0.5 

DM-derived distance (kpc) 8 

Most likely distance", d (kpc) 1 .5-3 



Note. — Numbers in parentheses represent the nominal 1 u TEMP02 
timing uncertainties on the last digits quoted. For gamma-ray parameters, 
the first uncer tainty is statistical and the second accounts for systematics (see 
Section r2.3.H . 

See discussion at the beginning of Section[3]and at the end of Section [3!2l 



